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Abstract 

The aim of this paper is to determine the locally densest horoball packing 
arrangements and their densities with respect to fully asymptotic tetrahedra 

3 

with at least one plane of symmetry in hyperbolic 3-space H extended with 
its absolute figure, where the ideal centers of horoballs give rise to vertices 
of a fully asymptotic tetrahedron. We allow horoballs of different types at 
the various vertices. Moreover, we generalize the notion of the simplicial 
density function in the extended hyperbolic space H , (n > 2), and prove 
that, in this sense, the well known Boroczky-Florian density upper bound 

3 

for "congruent horoball" packings ofH does not remain valid to the fully 
asymptotic tetrahedra. 

The density of this locally densest packing is « 0.874994, may be sur- 
prisingly larger than the Boroczky-Florian density upper bound « 0.853276 
but our local ball arrangement seems not to have extension to the whole hy- 
perbolic space. 



* Mathematics Subject Classification 2010: 52C17, 52C22, 52B15. 
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1 Basic notions 



1.1 Local density function 

We summarize the most important definitions and results about ball packings in 
H (n > 2). For more details and proofs, we refer to [fl], El, [HD> [f2 and ifToll . 
There are different notions of packing density. For later purposes, the local den- 
sity measure is the best suited one. In the n-dimensional (n > 2) hyperbolic 
space there are 3-types of spheres: sphere, horosphere and hypersphere. Now, 
we consider the horospheres and their bodies, the horoballs. A horoball packing 
Bh = {Bh} of H is an arrangement of non-overlapping horoballs Bh in H . The 
notion of local density of the usual ball packing can be extended for horoball pack- 
ings Bh of H . Let Bh G Bh, and P G H an arbitrary point. Then, p(P, Bh) is 
defined to be the length of the unique perpendicular from P to the horosphere Sh 
bounding Bh, where again p(P, Bh) is taken negative for P G Bh. The Dirichlet- 
Voronoi cell (shortly D-V cell) V(Bh) of Bh in Bh is defined to be the convex 
body 

V h = D(B h ,B h ) := {P G FT | p(P,B h ) < p(P,B' h ), \/B' h G B h }. (1.1) 

Since both, Bh and Vh, are of infinite volume, the usual concept of local density 
has to be modified. Let Q G <9H n denote the base point (ideal center at the infinity) 
of B h , and interpret Sh as a Euclidean [n — l)-space. Let B n _i(R) C Sh be an 
n — 1-ball with center C G Sh- Then, Q G <9H" and B n _ 1 (R) determine a convex 
cone C n (R) := cone(B < ^_ 1 (R)) G H" with apex Q consisting of all hyperbolic 
geodesies through B n _i(R) with limiting point Q. With these preparations, the 
local density 5 n (B h ,B h ) of B h to V h is defined by 

HR R , — vol(B h nC n {R)) 

5 n {B h ,B h ) := km , (1.2) 

R->ooVol{Vh H C n (R)) 

and this limes superior is independent of the choice of the center C of B n -i(R). 



1.2 Densest packings with horoballs of the same type 

We /zave to change the notion of the ,, congruent" horoballs in a horoball packing 
to the horoballs of the "same type " because the horroballs are congruent in the 
hyperbolic space H . Two horoballs in a horoball packing are in the "same type " 
if and only if the local densities of the horoballs to the corresponding cell (e.g. 
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D-V cell; or ideal simplex, later on) are equal. If we assume that the „horoballs 
belong to the same type", then by analytical continuation, the well known sim- 
plicial density function on H can be extended from n-balls of radius r to the 
case r = oo, too. Namely, in this case consider n + 1 horoballs B h which are 
mutually tangent. The convex hull of their base points at infinity will be a totally 
asymptotic or ideal regular simplex T reg G H 1 of finite volume. Hence, in this 
case it is legitimate to write 

d n (oo) = (n + l) y . (1.3) 

vol{T reg ) 

Then for a horoball packing there is an analogue of ball packing, namely (cf. 
Hi, Theorem 4) 

S n {B h , B h ) < d n (oo), \/B h e B h . (1.4) 

The upper bound d n (oo) (n = 2, 3) is attained for a regular horoball packing, that 
is, a packing by horoballs which are inscribed in the cells of a regular honeycomb 
of H . For dimensions n = 2, there is only one such packing. It belongs to the 
regular tessellation {oo, 3} . Its dual {3, oo} is the regular tessellation by ideal 
triangles all of whose vertices are surrounded by infinitely many triangles. This 
packing has in-circle density (^(oo) = - ~ 0.95493... 

In H 3 there is exactly one horoball packing whose Dirichlet-Voronoi cells give 
rise to a regular honeycomb described by the Schlafli symbol {6, 3, 3} . Its dual 
{3,3,6} consists of ideal regular simplices T reg with dihedral angle | building up 
a 6-cycle around each edge of the tessellation. 



1.3 Optimal packings by horoballs of different types in H 

In lIH we have determined the optimal horoball packing arrangements and their 
densities for all four fully asymptotic Coxeter tilings (Coxeter honeycombs) in 
H . Centers of horoballs are required to lie at vertices of the regular polyhedral 
cells constituting the tiling. We allow horoballs of different types at the var- 
ious vertices. We have proved that the known Boroczky-Florian density upper 
bound for "congruent horoball" packings of H remains valid for the class of 
fully asymptotic Coxeter tilings, even if packing conditions are relaxed by allow- 
ing horoballs of different types under prescribed symmetry groups. The conse- 
quences of this remarkable result are discussed for various Coxeter tilings (see 
[|9]|), but in this paper we consider only the tetrahedral Coxeter tilings. 
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a. b. 
Figure 1: Two optimal horoball arrangemets of (3, 3, 6) tiling. 



The Coxeter tiling {3, 3, 6} is a three-dimensional honeycomb with cells com- 
prised of fully asymptotic regular tetrahedra. We arbitrarily choose such a tetra- 
hedron T reg = E E 1 E 2 E 3 , and place the horoball centers at the ideal vertices 
E Q , . . . , E 3 . We vary the types of the horoballs so that they satisfy our constraints 
of non-overlap. The density of horoball packings to the above honeycomb can be 
computed by the new type of the simplicial density function: 



Definition 1.1 The density of a horoball packing in Coxeter honeycomb {3, 3, 6} 
is defined as 

J2t=i vol(Bi n T reg ) 



vol(T reg ) 



(1.5) 



Theorem 1.2 There are two distinct optimally dense horoball arrangements B\ m , 
(i = 1,2) for the tetrahedra! Coxeter tiling (3,3,6) with the same density: 5{B\ m ) 
0.85327609. 



Fig. 1. shows the two optimal horoball arrangements in Beltrami-Cayley-Klein 
model. 
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1.4 Generalization of the simplicial density function 

Definition 1.3 We consider an arbitrary fully asymptotic simplex T = E Ei 
E 2 E 3 . . . E n in the n-dimensional hyperbolic space H . Centers of horoballs 
are required to lie at vertices ofT. We allow horoballs (Bi, i = 1,2, ... ,n) of 
different types at the various vertices and require to form a packing, moreover we 
assume that 

card(Bi n [E io E h . . . £ in _J) < 1, ij ^ i, j E {0, 1, . . . ,n — 1}. 

(The hyperplane of points E io , E il , . . . , E in _ 1 is denoted by [E^E^ . . . E in _ 1 ] may 
touch the horoball B in .) The generalized simplicial density function for the above 
simplex and horoballs is defined as 



8(B) 



vol{T) 



For n = 3 the main problem is to find a fully asymptotic tetrahedron T opt e 

— 3 

H and horoballs Bi centered at the vertices Ei such that the density 5(B) (see 
Definition 1.3) of the corresponding horoball arrangement is maximal. In this 
case the horoball arrangement B is said to be locally optimal. 

In this paper, we study locally optimal horoball packings for fully asymptotic 

3 

tetrahedra in H , while allowing different types of horoballs to be centered at the 
vertices of the fully asymptotic tetrahedra. 

The general investigation of this problem seems to be very difficult, thus we 
assume for simplicity that the tetrahedron has at least one plane of symmetry, sim- 
ply called plane symmetric. So it will have only one free angle parameter. Then 
this ideal tetrahedron has 2 orthogonal symmetry planes and 2 symmetry lines of 
these planes, i.e. the symmetry group is of order 8, denoted by 2 *2 = 4m2 (by 
Conway and Hermann-Mauguin, respectively). 



2 Computations in projective model 

For H™ n > 2 we use the projective model in Lorentz space E 1,n of signature 
(1, n), i.e. E 1,n is the real vector space V n+1 equipped with the bilinear form of 
signature (1, n) 

(x, y) = -xy + xy + -- + xY (2.1) 
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where the non-zero vectors 

x = (x°, x 1 , . . . , x n ) G V" +1 and y = (y°, y 1 ,..., y n ) G V n+ \ 

are determined up to real factors and they represent points in V n (R). H n is rep- 
resented as the interior of the absolute quadratic form 

Q = {[ x ] G V n \( x, x) = 0} = dW 1 (2.2) 

in real projective space P n (V n+1 , V n+ \). All proper interior point x G H n are 
characterized by ( x, x) < 0. 

The points on the boundary dH n in V n represent the absolute points at infinity 
of H . Points y with ( y , y) > lie outside of H and are called outer points of 
H™. Let X([x]) G V n a point; [y] G V n is said to be conjugate to [x] relative to 
Q when ( x, y) = 0. The set of all points conjugate to X([x]) form a projective 
polar hyperplane 

pol(X):={[y]eV n \(x, y) = 0}. (2.3) 

Hence the bilinear form Q by (2.1) induces a bijection (linear polarity V n+1 —t 
V n+ i) from the points of V n onto its hyperplanes. 

Point X[x] and the hyperplane a[a] are called incident if the value of the 
linear form a on the vector x is equal to zero; i.e., xa = (x G V™ +1 \ 
{0}, a G V n +i \ {0}). Straight lines in V n are characterized by the 2-subspaces 
of V n+1 or (n — l)-spaces of V n+ i (see e.g. in |fTTfl). 

In this paper we set the sectional curvature of H , K = —k 2 , to be k = 1. The 
distance s of two proper points (x) and (y) is calculated by the formula: 

cosh s = ~\ x ' y) (2.4) 

v( x ' x )(y> y) 

The foot point Y(y) of the perpendicular, dropped from the point X(x) on the 
plane (u), has the following form: 

(x, u) 

y = x-- -u. (2.5) 

(u, u) 

2.1 On horospheres in hyperbolic space H 3 

Definition 2.1 A horosphere in the hyperbolic geometry is the surface orthogo- 
nal to the set of parallel lines, passing through the same point on the absolute 

3 

quadratic surface (simply absolute) ofH . 
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We represent H by the Beltrami-Cayley-Klein ball model. We introduce a 
projective coordinate system using vector basis bj (i = 0, 1, 2, 3) for V 3 where 
the coordinates of the center of the model is (1, 0, 0, 0). We pick up an arbitrary 
point at infinity to be £"3(1, 0, 0, 1) (see Fig. 3). 

We obtain the following equation for the horosphere in our Beltrami-Cayley- 

— 3 

Klein model of H' above: 

- 2sx°x° - 2x 3 x 3 + 2(s + l)(x°x 3 ) + (s - + x 2 x 2 ) = (2.6) 



Remark 2.2 1 . We get the equation of the horosphere in the Cartesian coor- 
dinate system (x : = y := z : = 

2. We will also use the equation of the horosphere with center E 2 (see Fig. 3): 

2(±x 2 + y 2 + \z 2 + &xz) A(fx - \z - i 



l-si (1- Sl ) 



2 



The length l(x) of a horocycle arc to a chord segment x is determined by the 
classical formula due to J. Bolyai: 

l{x)=ksinhy. (2.8) 
k 

The intrinsic geometry of the horosphere is Euclidean, therefore, the area A of 
a horospherical triangle will be computed by the formula of Heron. The volume 
of the horoball sectors can be calculated using another formula by J. Bolyai. If 
the area of a domain on the horoshere is A, the volume determined by A and the 
aggregate of axes drawn from A is equal to 

V = l -kA, (2.9) 



we assume that k — 1 in this paper. 
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2.2 The volume of the fully asymptotic tetrahedron in H 

A tetrahedron T in hyperbolic space is determined in general by its six dihedral 
angles i.e. the mutual angles between the four faces of this tetrahedron. In order 
to calculate the volume vol(T) of T; one can dissect T into six orthoschemes and 
their volumes can be calculated by the Lobachevsky formula. 

A plane orthoscheme is a right-angled triangle, whose area can be expressed 
by the well known defect formula. For three-dimensional spherical orthoschemes, 
Ludwig Schlafli, about 1850, was able to find the volume differential depending 
on differentials of the not fixed 3 dihedral angles. Already in 1836, Lobachevsky 
found a volume formula for three-dimensional hyperbolic orthoschemes O [J3j . 

The integration method for orthoschemes of dimension three was generalized 
by Bohm in 1962 to spaces of constant nonvanishing curvature of arbitrary 
dimension. 



Theorem 2.3 (N. I. Lobachevsky). The volume of a three-dimensional hyperbolic 
orthoscheme O C H is expressed with the dihedral angles a m , a\ 2) «23, (0 < 
®ij < f ) (Fig- 2) in the following form: 




A 



Figure 2: 



vol(O) = ~{£{a i + 9)- C{a m -9) + C(- + a 12 - 9) + 
+ £(| - a 12 -9)+ £(a 23 + 9)- C(a 23 -9)+2£(^-9)} 



where 9 6 [0, |) is defined by the following formula: 



tan(0) 




COS Ctoi cos a 23 
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and where £(x) :— — J log |2 sin t\dt denotes the Lobachev sky function. 

o 

We obtain the volume formula for the asymptotical orthoschems if A and A 3 
are at the absolute i.e. a := «oi = f — a i2 = «23 = 

voliO^) = ^C(a). (2.10) 

J. Milnor's volume formula of a fully asymptotic tetrahedron T(a ) (3) with dihe- 
dral angles a, (3, 7 is determined by its division into orthoschems (see f[T4l ): 

vol(T(a,(3))=£(a) + £(l3)+£(j), where a + (3 + 7 = tt. (2.11) 

As an easy consequence, we get the following 

Lemma 2.4 The volume formula vol(T(a)) of a fully asymptotic plane symmetric 
tetrahedron T(a) can be derived by the duplication law: 

vol(T(a)) = 2£(a) - £{2a) = -2£{a + |). (2.12) 

In (HI there are further results for the volumes of orthoschemes and simplices in 
higher dimensions. 



3 Horoball packings for asymptotic tetrahedra 

The aim of this section is to determine the optimal packing arrangement and its 

3 

densities for the fully asymptotic tetrahedra in H . 

We will make heavy use of the following Lemma (see also ifTTll ): 

Lemma 3.1 Let B\ and B2 denote two horoballs with ideal centers C\ and C2 
respectively. Take T\ and r 2 to be two congruent trihedra, with vertices C\ and 
Assume that these horoballs Bi(x) and B 2 (x) are tangent at point I(x) G C1C2 
and C1C2 is a common edge of the two trihedra T\ and r 2 . We define the point 
of contact 1(0) such that the following equality holds for the volumes of horoball 
sectors: 

V(0) := 2vol(B 1 (0) n n) = 2vol(B 2 (0) n r 2 ). 
If x denotes the hyperbolic distance between 1(0) and I(x), then the function 

V(x) := voliB^x) D n) + vol(B 2 (x) fl r 2 ) = Y^±( e ^ + e ~^) 
strictly increases as x — > ±00. 
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We arbitrarily choose a fully asymptotic tetrahedron T(oe) = E E 1 E 2 E 3 
with at least one plane of symmetry (see Section 1.4 and Fig. 3), and place the 
horoball centers at vertices E 0l . . . , E 3 . We vary the types of the horoballs so 
that they satisfy our constraints of non-overlap. The packing density is obtained 
by Definition 1.3. The dihedral angles of the above tetrahedron at the edges 
E Ei, E E 2 , E 1 E 3 , E 2 E 3 are denoted by a and it is clear that the dihedral 
angles at the remaining edges E E 3 and E X E 2 are n — 2a. 

We introduce a Euclidean projective coordinate system to the tetrahedron T(ce) 
by the following coordinates of the vertices: 



E := (1, 0, VT^, z), (-l<z< 1); E x = (l, 0, --); 

E 2 = (1,-^,0, -i); E 3 := (1,0,0,1). 

Eq lies on the symmetry plane x = of tetrahedron T(a) (see Fig. 3). In order to 
determine the optimal ball arrangement and its fully asymptotic tetrahedron first 
we consider the following situations: 

1 . In an optimally dense packing, the horoball packing must be locally stable, 
i.e. each ball is fixed by its neighboring horoballs or by the opposite face 
of the tetrahedron. Otherwise the density could be improved by blowing up 
at least one horoball until it touches a neighboring horoball or its opposite 
face of tetrahedron. 

2. We fix the parameter z (z £ (—1,1) \ 0) and blow up the horoballs Bi cen- 
tered at the vertices of T(a) so that the volumes vol(B{ fl T(a)), (i = 

0. 1.2.3) are equal until certain two neighbouring horoballs touch each 
other. In this situation we distinguish two cases for the varying vertex E : 

1. ) If < z < 1 then the horoballs B 3 — B and B\ — B 2 touch each other 
the other horoballs do not. 

ii.) If —1 < z < then each of the horoball pairs B — B 1} B — B 2 , Bi — 
£>3, B 2 — B 3 have exactly one common point, respectively, and B 3 — B 
and Bi — B 2 do not touch each other. 

If z = the tetrahedron is regular T(a) = T reg , (a = |) and each horoball 
touch all neigboring horoballs. This case is studied in details formerly in 
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3. The volume vol(T(a) is determined by Lemma 2.4, formula (2.12), and 
by the above machinery of hyperbolic geometry we obtain the following 
formulas: 



cos(2a) = — 
vol(T(a)) = -2C 



l + 2z 

—• ' 

1 

- arccos 



1<Z< 1), 

l + 2z\ 7T 

~ z-2 ) + 2 



(3.1) 



i/1: (z e(o,=% + 



a.) First, we define the tangent point 1(0) G E E 3 of horoballs -Bo(O) and B 3 (0) 
so that equalities hold for the volumes of all horoball sectors. 

The volume sum of horoball sectors -Bj(O) fl T(a), (i = 0, 1, 2, 3) in T(a) 
is 

4 

V(0) = J2( vol ( B i{0) n T{a))) = 4wo/(5 1 (0) n T(a)). 
i 

In the introduced coordinate system £?i(0) and B 2 (0) touch each other at 
(1, 0, 0, — |). Consider the point I(x) on the edge E Q E 3 where the horoballs 




Figure 3: 



Bi, (i — 0, 3) are tangent. Here x denotes the hyperbolic distance between 
1(0) and I(x). Let I(xi) £ E E 3 be such a point and parameter x 1 where 
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horoball B 3 (x±) touches -82(0) and -Bi(O). This horoball arrangement is 
denoted by B\{xi{z)). Function Vi(x) to B\ is defined by: 

V x {x) := ±V(0) + \v(0)e 2x + \v{0)e~ 2 \ x e (0,x 1 (z)}. 

By Lemma 3.1 it follows that function Vi(x) strictly increases as I(x) (x G 
(0, xi(z)}) moves away from 1(0) along E E 3 . That means that the density 
function for any z e (0, =f + (^f + « 0.64556191) 

^i(x) 



vol(T(a)) 

attains its maximum at the endpoint x\. Such a configuration does not occur 

if, e (-f + ^,l]. 

To study the former density function we have to compute the 

volume sum of the horoball sectors Vi(xi) = Ef (Sj fl T(a)) belonging to 
the parameter xi. Of course the volume Vi(xi) and the volume of the tetra- 
hedron T(a) depend on parameter z. First we calculate the six intersection 
points of the edges of T(cc) and horoballs Bi (i — 0, 1, 2, 3), using the above 
projective coordinate system and the equations of the horospheres derived 
from (2.7). These intersection points on the edge E^Ej is denoted by Ml 
(i<j), i,je {0,1,2,3}. 



Ml 



(3.2) 



/ y/l - z 2 2 + 2 \ , . V3 1\ 

( 1 '°'- 3 CT-2^)' ^ = d.-l-' '4)' 

^23= (l^' 'l)' M 12= (^0,0,-i), 

! _ / V3(z + 2) VT^ 5,-2 \ 
M ° 2 " I 1 ' 2(, + 5) ,d , + 5 '2(, + 5)J' 

! _ / V3> + 2) Q yr^ 5,-2 x 

01 ~ V ' 2(, + 5) ,d z + 5 '2(z + 5))- 

The volumes of horoball sectors (depending on parameter z) can be com- 
puted by Bolyai's formulas (2.8-9) and the volume of the tetrahedron T(a) 
is determined by formula (3.1). Finally, we obtain the density function 
8(Bi(xi(z))) = v ^t{1)) wn i cn depends only on parameter z e (0, +| + 
^]. By careful computer analysis of the above density function we get 
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that the function is convex, it attains its maximum at the upper endpoint 
of the interval (0, =| + ^3 ]• We have studied the above function with 
Maple using that the conditions of the Lobachevsky function are well known 
ll22l . The graph of 8(Bi(xi(z))) can be seen in Fig. 4. Note, that the case 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
0.0 0.1 0.2 0.3 0.4 0.5 0.6 



Figure 4: 

z = belongs to the optimal arrangement of Boroczky-Florian, its den- 
sity is ~ 0.85327609. The density at the upper endpoint of the interval 
(0, + ^] is larger than that, namely: 

6 { Bl { Xl {l5 + ~w))) ~ °- 86767481 ' ( a ~ i- 30899694 )- ( 3 - 3 ) 

b.) We consider the former horoball arrangement Bi(xi(z)) for any z G (0, ^|+ 

^p]. We expand the above "larger horoball" B 3 (xi) until it touches the op- 
posite face E E 1 E 2 of the tetrahedron while keeping the other horoballs 
tangent to it. The arising horoball arrangement are realizable for any z G 

v u ' 13 ' 13 J • 

Similarly to the case a.) we consider a point I(x) (x > X\) (on the edge 
E E 3 where the horoballs Bi, [i = 0, 3) are tangent at point I(x) G E E 3 . 
Let x — Xi denote the hyperbolic distance between I(xi) and I(x). Further- 
more, let I(x 2 ) G E E 3 be such a point where horoball B 3 (x 2 ) of parameter 
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x 2 touches the face E E 1 E 2 . The "largest horoball" determines the config- 
uration of all other horoballs and this horoball arrangement is denoted by 
B 2 (x 2 (z)). Function V 2 (x) is defined as follows: 



The following Lemma is obtained by examining the above function: 

Lemma 3.2 The maxima of function V 2 (x) are realized to parameters x\ or 
x 2 . 

By the above Lemma 3 .2 it follows that it is sufficient to consider the volume 
function V 2 (x) and the density function 5(B 2 (x(z))) at the parameters x\ 
and x 2 . The case x\ was taken in a.) thus we have to consider only the case 

x = x 2 . 

The horoball B 3 (x 2 ) has to touch the side face E E 1 E 2 of fully asymptotic 
tetrahedron T(a). Thus, it is passing trough the foot point E' 3 perpendicu- 
larly dropped from the point E 3 on the plane E EiE 2 , 



To examine the density function S(B 2 (x 2 (z))) we have to compute the sum 
of the volumes of the horoball sectors V 2 (x 2 ) = £f(-B, n T(a)) to the 
parameter x 2 . The volume V 2 (x 2 ) and the volume of the tetrahedron T(a) 
depend on parameter z. Similarly to the case a.) we have to determine 
the seven intersection points of the edges of T(a) and horoballs Bi (i = 
0, 1, 2, 3). The intersection points on the edges E Ei, E E 2 , E E 3 , E 3 E 2 
and E 3 E 1 are denoted by Mf- (i < j), i,j e {0, 1,2,3}. On the edge E X E 2 
there are no touching point but there are two point of intersections M\ 2 . 




x e [x 1 ,x 2 \. 





(3.4) 
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1) VI - z 2 (z + 2) 3z 2 + 2z-2^ 



Az 2 — z — 6 



4z 2 — z — 6 



2V3(^ 2 



l) VT^(^ + 2) 



4z 2 — 2 — 6 J 
3z 2 + 2z-2 



Az 2 



6 



Az 2 



6 



4z 2 



6 



The volumes of horoball sectors (depending on parameter z) can be com- 
puted by Bolyai's formulas (2.8-9) and the volume of the tetrahedron T{a) 
is determined by formula (3.1). We get the density function 8(B 2 (x2(z))) = 
vd{T(a)) wn i cn depend only on parameter z G (0, -=| + By careful in- 
vestigation (see 11221 ) of the above density function we get that the function 
attains its maximum at the point z = =j + The graph of 8(B 2 (x2{z))) 
is shown in Fig. 5. The density at the point ^ 



^j^- is larger than the 



0.865- 



0.86- 



0.855- 



0.85- 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
0.0 0.1 0.2 0.3 0.4 0.5 0.6 



Figure 5: 

Boroczky-Florian upper density, namely 

5(b 2 (x 2 (^ + ^p))) » 0.86767481, (a « 1.30899694). (3.5) 

The horoball density of of Boroczky-Florian belongs to the parameter 2 = 
0. To the parameter z = ^| + ^ we get the same horoball density as in 
case a.). 

i/2: (ze + M 1]) 
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Figure 6: 

a.) ze , where z 3 := ^y3 v / T78cos(|arctan(gg T v / 3v / 229))- 

f| « 0.9061774494. 

In the former case i/l/b we have already considered the ball arrangement 
Bi(xi(z)) for any z G (0, ^ + ]• There we have expanded the "larger 
horoball" B 3 (xi) until it touches the opposite face E EiE 2 of the tetra- 
hedron while keeping other horoballs tangent to it. But if z G (-jf + 
^^,z 3 ] then we cannot apply this procedure. Namely, if the horoball 
B 3 (x 2 ) touches the opposite side face, then B 3 (x 2 ) does not touch 2?i(0) 
and -82(0). Then we expand the horoball B 2 until it touches B 3 (x 2 ) or the 
opposite side face. This means that the "larger horoball" .83(0:2) touches 
the opposite face E EiE 2 at point E' 3 and the horoballs B 2 , B touch B 3 . 
Then Bi touches only the horoball B 2 . This configuration occurs until the 
horoball B 2 is not tangent to its opposite face E 3 E 1 E . 

This arrangement exactly exists if z G (-jf + ^j-, 23] • The horoball B 2 
touches the opposite side face E E 1 E 3 of T(a) at the point E 2 : 

= / 5^ + 7 ^3(2;? + 1) 3y^? 2z + l \ 
2 V2(^ + 2)' 2(^ + 2) ' 2(0 + 2) ' 2(z + 2)A 

Similarly to the case a.) the "largest horoball" determines the configuration 
of all other horoballs and this horoball arrangement is denoted by B3 (x 3 (z) ) . 
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The function V 3 (x 3 (z)) and related density function 

Vs(x 3 (z)) 



vol(T(a)) 



can be analysed in the same way as in i/l/a and i/l/b and we obtain that the 
convex density function attains its maximum at the above z 3 . 

The graph of 5(B 3 (x 3 , z)),z G + ^$ )> ^3] i s shown in the Fig. 6. While 
Fig. 7 shows the larger horoballs B 3 , B 2 E B 3 (x 3 (z 3 )) with tetrahedron 
T(a) in the Beltrami-Cayley-Klein model. The density at the point z = z 3 



!MW5iBBBBBiSSii£8!l*2S£ 
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Figure 7: The optimal horoball packing arrangement related to a fully asymptotic 
tetrahedron 

is larger than the above maxima 

5(B 3 {x 3 (z 3 ))) ps 0.87499429, (a ps 1.44340117). (3.7) 

b.) z e (z 3 , 1) 

We consider such ball arrangements B^x^z)) where the "larger horoballs" 
B 3 and B 2 touch their opposite faces E EiE 2 and EqE\E 3 , respectively. 
The horoballs B , B 3 and B±, B 2 pairwise touch each other. This horoball 
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packing can be studied similarly to the above cases and we can analyse the 
density function 8(B±(x±(z))) (see Fig. 8), not detailed further. We obtain, 
that the density function is stricly decreasing in the interval (z 3 , 1). Thus, 
there is no horoball arrangement with larger density than B 3 (x 3 (z 3 )) with 
density « 0.87499429 (see (3.7)). 




Figure 8: 



Analogously to (i) we have investigated the cases (ii) (—1 < z < 0), not detailed 
more, and we have obtained the following 

Theorem 3.3 The densest horoball arrangement, related to the generalized sim- 
plicial density function (Definition 1.3) belongs to the horoball arrangement B 3 (x 3 (z 3 )) 
(see (3.7) and Fig. 7) with density 

8(B 3 (x 3 (z 3 ))) « 0.87499429. 

Remark 3.4 This packing seems not to have global extension to the entire hyper- 
bolic space H , so that the same density occurs in each asymptotic tetrahedron. 
We plan to investigate this problem with E. Molndr and I. Prok by /fi2~l/. 

Optimal sphere packings in other homogeneous Thurston geometries repre- 
sent a class of open mathematical problems. For these non-Euclidean geometries 
only very few results are known [|20l . f|2T|. Detailed studies are the objective of 
ongoing research. 
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